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The Sunyaev-Zel'dovich (SZ) effect of a cluster of galaxies is usually measured after 
background radio sources are removed from the cluster field. Gravitational lensing by 
qq ■ the cluster potential leads to a systematic deficit in the residual intensity of unresolved 

sources behind the cluster core relative to a control field far from the cluster center. 
As a result, the measured decrement in the Rayleigh-Jeans temperature of the cosmic 
microwave background is overestimated. We calculate the associated systematic bias 
| which is inevitably introduced into measurements of the Hubble constant using the SZ 

effect. For the cluster A2218, we find that observations at 15 GHz with a beam radius 
Qh' of 0'.4 and a source removal threshold of 100^ Jy underestimate the Hubble constant 

Q ■ by 6-10%. If the profile of the gas pressure declines more steeply with radius than 

» ■ that of the dark matter density, then the ratio of lensing to SZ decrements increases 

Ci ■ towards the outer part of the cluster. 

> 

X 
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1. Introduction 

The Sunyaev-Zel'dovich (SZ) effect describes the distortion introduced to the Cosmic 
Microwave Background (CMB) spectrum due to its Compton scattering off free electrons, which 
are either hot (the thermal effect) or possess a bulk peculiar velocity (the kinematic effect; see 
reviews of both effects in Sunyaev & Zel'dovich 1980 and Rephaeli 1995). The thermal SZ effect 
provides an important diagnostic of the hot gas in clusters of galaxies, and by now has been 
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measured in a number of systems (see Table 1 in Rephaeli 1995). The kinematic effect has an 
amplitude which is typically an order-of magnitude smaller and has not yet been definitively 
detected (see Rephaeli & Lahav 1991 and Haehnelt & Tegmark 1995, regarding prospects for a 
future detection). 

It has long been realized that a measurement of the thermal SZ effect, combined with X-ray 
observations, can be used to estimate the distance to the cluster and hence the Hubble constant, 
Hq, under the assumption that the cluster is spherical (Cavaliere, Danese, & De Zotti 1977, Gunn 
1978, Silk & White 1978, Birkinshaw 1979). The inferred value of the Hubble constant is inversely 
proportional to the square of the SZ temperature decrement. This approach had led to values of 
the Hubble constant which are typically on the low side of the range inferred from other methods 
(see, e.g., Table 2 in Rephaeli 1995). An often cited systematic effect that could account for 
this bias is elongation of the selected clusters along the line-of-sight. In this Letter, we explore a 
different effect which leads to a systematic bias towards Iow-Hq values even if these clusters are 
perfectly spherical. The effect results from gravitational lensing by the cluster potentials. 

Measurements of the decrement in the Rayleigh- Jeans (RJ) temperature of the microwave 
background due to the thermal SZ effect are routinely accompanied by the removal of background 
radio sources down to some flux threshold (see, e.g., Birkinshaw, Hughes, & Arnaud 1991). In 
this process, it is implicitly assumed that the flux threshold for the removal of sources behind the 
cluster core is the same as in a control field far from the cluster center. However, this assumption 
is not strictly true due to the inevitable magnification bias which is introduced by the gravitational 
lensing effect of the cluster potential. In reality, the cluster acts as a lens which magnifies and 
thus resolves sources that are otherwise below the detection threshold. The residual intensity of 
unresolved sources is therefore systematically lower behind the cluster core, as compared to that 
in the control field. Lensing artificially increases the flux deficit behind the cluster core and thus 
leads to a systematic underestimate of the Hubble constant. 

In this Letter we calculate the effect of lensing on SZ measurements of the Hubble constant. 
Our discussion on lensing follows closely the approach developed in an earlier paper (Refregier &: 
Loeb 1996, hereafter RL) which focused on lensing of the X-ray background by galaxy clusters; 
the interested reader should consult this earlier paper for more details. Here, we describe our 
models for the background population of radio sources and for the cluster potential in §||. We then 
show in §H] how the lensing effect leads to a systematic decrement in the intensity of unresolved 
sources. In §|I|, we present numerical results for different values of our model parameters and for 
the specific example of A2218. Finally, §|| summarizes the main conclusions of this work. 

2. Model 

We model the gravitational potential of the cluster as a Singular Isothermal Sphere (SIS) (e.g. 
Schneider et al. 1992). This model provides a good first-order approximation to the projected 
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mass distribution of known cluster lenses (Tyson & Fischer 1995, Narayan & Bartelmann 1996, 
Squires et al. 1996a,b). The SIS potential causes background sources to appear brighter but 
diluted on the sky by the magnification factor 

(1) 

where 6 is the angle between the image of the source and the lens center, and a is the Einstein 
angle. For a SIS with a line-of-sight velocity dispersion a v , a = 47r(a"^/c 2 )D/ s /L> os , where Di s and 
D os are the lens-source and the observer-source angular diameter distances, respectively. 

In general, the Einstein angle a depends on the redshifts of the lens, Z[, and of the source, 
z s . However, the dependence on the source redshift is weak if z s ;> 3zi (cf. Fig. 1 in RL). Most 
measurements of the SZ effect are performed with nearby clusters (zi <^ 0.2), while sub-mJy radio 
sources have median redshifts in the range 0.5 <^ z s ^ 0.75 (Windhorst et al. 1993). We therefore 
take a to be independent of z s and simply consider the two-dimensional distribution of radio 
sources on the sky. 

We model the flux distribution of background radio sources according to the observed 
number-flux relation at 4.86 GHz (Windhorst et al. 1993). The squares in Figure [l] show the mean 
values of the observed differential counts, dn/dS, normalized by S~ 2 . The dotted line illustrates 
the number counts limits inferred by Fomalont et al. (1991) from a fluctuation analysis of the radio 
source background (see also Windhorst et al. 1993). The observed counts extend from ~ 10~ 5 
to 10 Jy. We model these counts by six broken power laws of the form, (dn/dS)\ s = r]jS~ Xj 
for Sj-ij > S > Sjj + i, with j = 1, ... ,6. The six slope parameters are Ai = 2.50, A2 = 2.64, 
A3 = 2.34, A4 = 1.67, A5 = 2.19, and A6 = 1.70. The power-law break points are S12 = 10, 
S 23 = 6.0 x 10 _1 , S34 = 2.6 x 10~ 2 , S45 = 8.0 x 10~ 4 , and S 56 = 1.2 x 10~ 6 Jy. The normalization 
is set by imposing continuity and using rj\ = 54.0 Jy Al_1 sr _1 . The model relation is shown as 
the solid line in figure |]. The counts were extended below 10 ^Jy based on the extrapolation 
suggested by Windhorst et al. (1993). For a typical SZ threshold Sa S> 10/iJy, the results are not 
sensitive to the parameters of this extrapolation. Our total 4.86 GHz intensity from radio sources 
is i to4 (4.86GHz) = / °° S(dn/dS)dS = 3.7 x 10 3 Jy sr" 1 or equivalently 5.1 mK. 

At frequencies above 4.86 GHz, we approximate the mean spectrum of the sources by a power 
law, S u oc v~~< . Based on the observed fluxes between 1.41-8.44 GHz, Windhorst et al. (1993) 
estimate the median spectral index of sources with S4.86 ~ 0.1 mJy to be 7 = 0.35 ± 0.15. Radio 
sources with S^&% ^ 1 mJy have a median angular size ^ 2" (cf. Fig. 2 in Windhorst et al. 

1993) . In the FIRST radio survey, White et al. (1997) found that only ~ 20% of all sources with 
S^l^GHz) ^ 1 mJy have a major axis ^ 5". The angular extent of sub-mJy radio sources is thus 
much smaller than the Einstein angle of clusters with observed arcs (a ~ 30"; see Le Fevre et al. 

1994) , as well as the beam size used in SZ measurements (^ 1'; cf. Rephaeli 1995). Thus, we can 
safely ignore the finite extent of the radio sources and treat them as if they were pointlike. 
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3. The Lensing Effect 



In a region of the sky where the magnification factor is /x, the apparent differential count of 
sources obtains the value (dn/dS)\ s = [i~ 2 (dn/dS) ^ , where hat denotes unlensed quantities. In 

particular, for a power law differential count relation, (dh/dS) , oc S~ x , the observed differential 
count is, (dn/dS)\ s oc n^~ 2 S~ x . The differential count therefore increases (decreases) as [i 
increases if A is above (below) the critical slope A cr it = 2. When A = A cr it, lensing has no effect 
on the apparent differential count. Interestingly, figure |] shows that the actual radio count slope 
oscillates around A cr i t for fluxes S4.86 ^ Uy- 

In measurements of the SZ effect, discrete sources are typically removed down to a given 
detection flux threshold, Sd- The mean residual intensity i(< Sd) due to the superposition of all 
undetected discrete sources with fluxes below Sd is then assumed to be equal to its sky-averaged 
value 

»(< S d ) = dSS^r 

Jo dS 

However, the magnification due to lensing lowers the unresolved intensity systematically relative 
to its sky-averaged value and changes it to, 



(2) 

s 



l (<s d )= j Sd dss^^ 



= i(< S d /n), (3) 



where \i = fi(9) is given by equation (|T]). Lensing conserves the total intensity of the radio 
source background and merely reduces the effective flux threshold for resolving sources by a 
factor fi. The intensity offset due to lensing is then, Aii ens = i(< Sd) — i(< Sd)- In the RJ 
regime, this can be expressed more conveniently in terms of the brightness temperature difference, 
ATJcns = (c 2 /2v 2 k-Q) Azi ens , where fee is Boltzmann's constant. For \i > 1 (i.e., 6 > a/2), the 
unresolved intensity is decreased, implying a negative ATi ens , and so the SZ decrement in the 
RJ regime, ATsz, is overestimated due to lensing. Note that for 6 S> a, equation (|l]) yields 
fj, m 1 + a/9 and ATi ens cx O^ 1 . 

The effect of lensing on estimates of Hq can be easily found from the scaling, Hq cx (ATsz)~ 2 , 
where ATsz is the temperature offset produced by the SZ effect. The small systematic correction 
(Ai^o)icns = Ho(tiue) — Hq (observed), which must be incorporated in order to compensate for the 
lensing effect is, to leading order, 

(Aff )i cns ^ ^ ATi ens 

Hq ATsz 

For 6 ^ a/2 and the RJ spectral regime, both ATi ens and ATsz are negative and so (AiTo)iens 
is positive. The lensing correction will then tend to increase the estimated value of the Hubble 
constant. 
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4. Results 

Figure ^ shows ATi ens as a function of angular separation from the cluster center, 9, for several 
values of the detection threshold Sd- The values for ATi ens and Sd correspond to a frequency of 
4.86 GHz. The dependence of A7] ens on the Einstein angle a and on frequency v were conveniently 
factored out. The Einstein angles for clusters with observed optical arcs are in the range of 10-50" 
(Le Fevre et al. 1994). 

The lensing decrement, ATi ens , shows a sharp peak near the Einstein angle. For 9 a, 
ATi ens is first weakened and then enhanced as Sd varies from 10 -2 to 10 -5 Jy. This is due to the 
fact that the count slope A crosses the critical value A cr jt = 2 around Sd ~ 10~ 3 Jy (see Fig. ||). 
The enhancement in ATi ens as Sd decreases below 10~ 3 Jy occurs in spite of the decrease in the 
unlensed intensity i(< Sd) there. At these fluxes, the removal of fainter radio sources paradoxically 
makes the lensing decrement more pronounced. Note that because of the large shot noise in the 
source counts (with an rms of at/i rj 0.5 in a 1 arcmin 2 cell for Sd(4.86GHz) = 10 -3 Jy), ATi cns 
will not necessarily be realized in each individual cluster. The lensing induced decrement should 
be regarded as a systematic effect that must be corrected for statistically, when a large sample 
of clusters is considered. For observations with a large field of view, the lensing signature might 
appear in the outer part of each individual cluster. 

As a specific example we consider A2218, an Abell richness class 4 cluster at a redshift 
Z\ = 0.175, which shows several optical arcs (Pello et al. 1992, Le Borgne et al. 1992). Arc no. 359 
in Pello et al. (1992) is separated by 20". 8 from the central cD galaxy and has a measured redshift 
of 0.702, close to the probable median redshift of sub-mJy sources (z s ~ 0.5-0.75; cf. Windhorst 
et al. 1993). We therefore model the cluster potential as a SIS with an Einstein angle of a = 20". 8 
for our radio sources (see also Miralda-Escude & Babul 1995). Interferometric imaging of the SZ 
effect in this cluster was performed by Jones et al. (1993) at 15 GHz, after the removal of point 
sources with fluxes above Sd(15GHz) ~ 1 mJy The restoring beam for their short baseline image 
had a FWHM of 129" x 120". The observed angular dependence of ATsz was fitted by a /3-model, 

ATsz(#) = ATo(l + ^j . (5) 

Acceptable x 2 values were obtained for different sets of parameters ranging from (3 ~ 0.6, 9 C 0'.9, 
and AT « 1.1 mK, to /3 « 1.5, 9 C « 2'.0, and AT = 0.6 mK. 

Figure ^ shows the expected ratio ATi ens /ATsz for A2218 at 15 GHz, assuming a source 
spectral index of 7 = 0.35. The ratio is shown for two values of Sd(4.86GHz) and for the two 
extreme sets of fit parameters for ATsz(#)- The sharp peak at 9 = 20". 8 reflects the enhancement 
in ATi ens around the Einstein radius (see Fig. |2|). For the f3 = 1.5 model, ATi ens /ATgz diverges at 
9 £ 2'. If the mass distribution follows the SIS profile, AT lens oc 9' 1 at 9 » a. Since AT SZ oc 9 1 - 313 
for > 9 C (cf. Eq. §), the ratio (AT lens /AT sz ) oc 6» 3 ^" 2 diverges at large radii if fi > 2/3. The 
values of (3 derived from X-ray observations of clusters have a large scatter around a mean value 
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of ~ 0.65 (Sarazinl988, Jones & Forman 1984, Bahcall & Lubin 1994). Weak lensing studies in the 
optical band could be used in conjunction with X-ray observations to predict the relative radial 
behavior of the lensing and SZ decrements in each individual cluster. 

It is convenient to average the temperature offset over a circular "top hat" beam of radius 
Of, centered on the cluster center, (AT (8b)) = 28^ 2 Jq 6 8AT(8)d8. For the above model of 
A2218 with Sd(4.86GHz) = 10~ 4 Jy, the 15 GHz mean temperature offsets due to lensing 
are (ATi ens ) « —28,-14, and —1.6 [xK, for 8b = 0.4,1, and 60 arcmin, respectively. The 
corresponding decrement ratios are (ATi ens )/(ATsz) ~ 0.05,0.03, 0.18, for the j3 = 1.5 fit, and 
0.03,0.02,0.002 for the f3 = 0.6 fit. The fractional correction to the Hubble constant (Eq. |j§) is 
then AHq/Hq ~ 6-10%, 4-6%, and 0.4-40% for 8b = 0.4, 1, and 60 arcmin, respectively, where the 
ranges reflect the ambiguity in the fit parameters of AT$z(8). 



5. Conclusions 

We have shown that gravitational lensing of unresolved radio sources leads to a systematic 
overestimate of the SZ temperature decrement at angles 8 > a/2. The amplitude of the lensing 
effect peaks close to the Einstein angle of the cluster, a ~ 30" (cf. Fig. 2). While ATsz is 
independent of frequency in the RJ regime, the lensing decrement ATi ens oc z;~ 2 ~ 7 (with 7 ~ 0.35) 
is significant only at frequencies v <^ 30 GHz. In clusters where the radial profile of the gas pressure 
is steeper than that of the dark matter density (e.g., due to a gradient in the gas temperature), 
the ratio of the lensing to the SZ decrement increases at large projected radii. For observations 
of A2218 at 15 GHz with a source removal threshold of «Sd(4.86GHz) = 10 -4 Jy, Hq could be 
overestimated by ~ 0.4-40%, for a beam radius in the range of 0.4-60 arcminutes (cf. Fig. 3). 
The importance of the lensing effect will be enhanced in future observations (including attempts 
to detect the kinematic SZ effect) with greater sensitivity, higher angular resolution, and fainter 
source removal threshold. 

Lensing should also affect the power spectrum of microwave background anisotropies on 
~ 1' scales behind the cluster. These anisotropies are expected to originate primarily from the 
Ostriker-Vishniac effect (Hu & White 1996) and the cumulative SZ effect of other background 
clusters (Colafrancesco et al. 1994, Rephaeli 1995). Future SZ experiments might be contaminated 
by noise from these fluctuations (AT/T <^ 10 -6 ), especially in the outer parts of clusters. However, 
since these diffuse fluctuations will not be removed, lensing will conserve their net intensity and 
will not systematically offset the SZ decrement as it does in the case of discrete sources. 

We thank D. Helfand for useful comments on the manuscript. This work was supported in 
part by the NASA grants NAG5-3085 (for AL) and NAGW2507 (for AR). 
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Fig. 1. — Number-flux relation for radio sources at 4.86 GHz. The counts were normalized to S~ 2 , 
the relation which remains invariant under lensing. The approximate mean counts summarized by 
Windhorst et al. (1993) are shown as squares. The dotted line corresponds to the limits inferred 
from a fluctuation analysis of the unresolved background (Fomalont et al. 1991). The solid line 
shows our model with its six power-law components. 
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Fig. 2. — Temperature offset ATi ens induced by gravitational lensing as a function of angular 
separation from the cluster center. The offset is shown for several values of the source detection 
threshold S^. Values for A7] cns and correspond to a frequency v = 4.86 GHz. The dependence 
of ATi ens on the Einstein angle of the cluster, a, and on frequency, ^, were factored out. The 
parameter 7 is the mean spectral index of the radio sources, and ^4.86 = ^/(4.86GHz). 
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Fig. 3. — Lensing effect in A2218. The ratio of the temperature decrement induced by lensing at 15 
GHz, ATJens, to that induced by the SZ effect, ATsz, is shown for two values of the flux detection 
threshold S , ^(4.86GHz). The decrement ratio is evaluated for the two extreme fits obtained by Jones 
et al. (1993) to the observed radial dependence of ATsz- The curves which converge (diverge) at 
large radii correspond to the (3 = 0.6 ((3 = 1.5) fit. 



